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Abstract — Both the first-order signal statistics (e.g. the outage 
probability) and the second-order signal statistics (e.g. the av- 
erage level crossing rate, LCR, and the average fade duration, 
AFD) are important design criteria and performance measures 
for the wireless communication systems, including the equal 
gain combining (EGC) systems in presence of the cochannel 
interference (CCI). Although the analytical expressions for the 
outage probability of the coherent EGC systems exposed to CCI 
and various fading channels are already known, the respective 
ones for the average LCR and the AFD are not available in 
the literature. This paper presents such analytical expressions 
for the Rayleigh fading channel, which are obtained by utilizing 
a novel analytical approach that does not require the explicit 
expression for the joint PDF of the instantaneous output signal- 
to-interference ratio (SIR) and its time derivative. Applying 
the characteristic function method and the Beaulieu series, we 
determined the average LCR and the AFD at the output of 
an interference-limited EGC system with an arbitrary diversity 
order and an arbitrary number of cochannel interferers in forms 
of an infinite integral and an infinite series. For the dual diversity 
case, the respective expressions are derived in closed forms in 
terms of the gamma and the beta functions. 

Index Terms — Level crossing rate, average fade duration, 
cochannel interference, equal gain combining, Beaulieu series, 
Rayleigh fading. 



I. Introduction 

EQUAL gain combining (EGC) is an important diversity 
technique that is often used to mitigate fading in various 
wireless communications systems [1]. The EGC has several 
practical advantages over other diversity techniques, because 
it has close to optimal performance and yet is simple to imple- 
ment. The outage probability (OP) is the primary performances 
measure for all diversity systems, particularly those exposed 
to cochannel interference (CCI) such as the cellular mobile 
systems. When the CCI is the predominant noise source at 
the receiver, its OP represents the first-order statistical property 
of the output signal-to-interference ratio (SIR). The OP of the 
interference-limited EGC systems was studied in [2]-[3] and 
references therein. Apart from the OP, some aspects in the 
design and the analysis of wireless communication systems 
must also consider the signal correlation properties, thus 
necessitating the determination of its second-order statistical 
properties: the average level crossing rate (LCR) and the av- 
erage fade duration (AFD). They are used for proper selection 
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of adaptive symbol rates, interleaver depth, packet length and 
time slot duration in various wireless communication systems. 
While these statistics have already been determined for the 
signal envelope at the output of EGC systems exposed to 
various fading channels and thermal noise [4] -[6], the SIR 
statistics of the EGC systems subject to CCI have not yet 
been derived analytically to the best of author's knowledge. 
The average LCR and the AFD of the SIR at the output 
of selection combining (SC) and maximal-ratio combining 
(MRC) systems exposed to CCI and various fading (Rayleigh, 
Rice and Nakagami) channels have been reported only re- 
cently in [7]-[8], but these works do not consider the EGC 
systems. This paper focuses specifically on an interference- 
limited coherent EGC system with an arbitrary diversity order 
and an arbitrary number of cochannel interferers, and derives 
analytical solutions for the average LCR and the AFD of the 
output SIR in Rayleigh fading channels. 

Section II presents the coherent EGC system model and 
the channel model, including the two feasible scenarios for 
interference combining. Section III presents the analysis that 
yields to the analytical solutions of the OP, the average LCR 
and the AFD in forms of an infinite integral and an infinite 
series. The Section IV compares the computational burden 
between these two solutions and provides several numerical 
examples that illustrate the behaviors of the first-order and the 
second-order signal statistics. Section V summarizes the main 
results and concludes the paper. 

II. System and Channel Models 

We consider a coherent EGC communication receiver with 
M diversity branches. It is exposed to the transmissions 
of a single desired and N interference users, whose signal 
replicas in each diversity branch are received over independent 
identically distributed (IID) Rayleigh flat fading channels. 

In each diversity branch k (1 < fc < M), the desired 
signal is assumed to have an average power ils, while all 
interference signals have an equal average power fti. Thus, the 
channel gains in each branch can be represented as equivalent 
complex zero-mean Gaussian random variables (RVs) Wi^k', 
more particulary, Wo^k ~ -^^o.fc e^""''' with variance fts repre- 
sents the desired signal in branch k, while Wi,k — Yi^k e^^^-'' 
(1 < i < N) with variance f^i represents i-th interference 
signal in branch k. The phases of the desired signals Oi^^k 
and the interference signal 9i,k follow the uniform probability 
distribution function (PDF) over [0, 2tt), while the respective 
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envelopes Xo,k = y/\Wo,k\'^ and Fj,fe = -y/|Wi,feP follow the 
Rayleigh PDF. 

Due to the transmitter/receiver mobility and their relative 
velocity, the fading channel introduces time correlation of the 
real and imaginary parts of Wo,k (i-c, in-phase and quadrature 
components of the desired signal) with maximum Doppler 
frequency shift fmO in their power spectra. Additionally, the 
real and imaginary parts of the channel gains Wi^k of each 
interfering signal i are also assumed be time correlated with 
an identical maximum Doppler frequency shift /mi- 

In EGC systems, the desired signal replicas in each of the M 
branches are co-phased, equally weighted, and then coherently 
added to give the resultant desired output signal. For the 
interference combining, there are two possible scenarios: the 
signal replicas originating from any interferer can combine 
either incoherently [2, Section III] or coherently [3]. 

A. Incoherent Interference Combining 

If the interference signals are combined incoherently at the 
EGC output, the instantaneous SIR Z\ is determined as [2, 
Eq. (9)], 



M 



(1) 



where the powers (i.e. the squared envelopes) of all inter- 
ference signals in all diversity branches are added together. 
Thus, any single element in the denominator of (1), 1^^^, is 
a chi-squared RV with 2 degrees of freedom, so the entire 
denominator, given by 

N M 

^i' = EE^4' (2) 

i=l fe=l 

follows the chi-squared PDF with 2MN degrees of freedom. 
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where r( ) is the gamma function, defined by T{a) = 
j;^t^-^e-*dt [14]. 

B. Coherent Interference Combining 

If the interference signals are combined coherently at the 
EGC output, the instantaneous SIR Z2 is determined as [3, 
Eq. (2)], 
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where it is assumed that the roll-off factor of the equivalent 
baseband communication system is zero. In (4), the complex 
interference signals from all branches are first added together 
and then squared. Thus, J2iLi Wi,k is a complex Gaussian RV 
with zero mean and variance MOj, while its squared envelope 
(i.e. its power) is a chi-squared RV with 2 degrees of freedom. 
Thus, the denominator in (4), given by 
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follows the chi-squared PDF with 2N degrees of freedom. 
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The PDF of the numerators in (1) and (4), which are the square 
of the desired output signal envelope 
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(7) 
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is not known in closed form, except for M = 2. Thus, we 
revert to using the characteristic function (CF) method to 
arrive at the desired results. 

III. Average LCR AND AFD 

A. Definitions 

We first concentrate on the RV defined as the ratio of 
the envelopes of the desired signal X and the equivalent 
interference signal Yi - for incoherent interference combining, 
and Y2 - for coherent interference combining. 



Gi = 



and 



Go = \/ Zn = 



(8) 



(9) 



and denoted as the instantaneous envelopes ratio. We will 
first establish the average LCR of the envelopes ratio G and 
then readily obtain the average LCR and AFD of the SIR 
Z based on (8)-(9). The average LCR of the envelopes ratio 
G at threshold g is defined as the rate at which the fading 
process crosses level g in the negative direction [1]. It is 
mathematically defined by the Rice's formula [1, Eq. (2.106)] 



Ng{9) 



yfociy, 9)dg, 



(10) 



where G denotes the time derivative of G, and fQQ{g,g) is 
the joint PDF of G and G. The AFD is defined as the average 
time that the envelopes ratio G remains below the level g after 
crossing that level in the downward direction, and is defined 
by 

Fg{9) 



Tg{9) = 



Ng{9) ' 



(11) 



where Fg{ ) denotes the cumulative distribution function 
(CDF) of G. Considering (8)-(9), we introduce 9 = -Jz into 
(10) and (11), and determine the average LCR and the AFD 
for the SIR Z at threshold z as Nz{z) = Ng{^/z) and 
Tz{z) = Tg(v^), respectively. 

B. Characteristic Functions 

The desired signal at the EGC output consists of M 
Rayleigh RVs Xo.fc, each having as average power O = Os- 
Although the PDF of X is not known in closed form (except 
for M = 2), it is still possible to determine its CF in terms of 
the CFs of the constituent Xo,feS. 
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For this purpose, we define some general (Nakagami-like) 
RV U with the PDF given by (u > 0), 



1 \ " 2u 



" 0„,2a-l 



T{a) 

whose CF is given by [2] 



■ exp 



u 



(12) 



where coq = {2'k/T), T is a parameter governing the sampling 
rate in the frequency domain and controls the accuracy of 
the result, and Ai is an error term that tends to zero for 
large T. We assume T is large enough to omit this error 
term. Introducing (17) over (14) and changing the orders of 
summation and integration, we obtain 



$(a;,0,a)1l^J^°°/(«)exp(ia;«)dn=iF,(a;i;-^^fi) Fa{9) = \- £ 2lu.{^^{n.,)^*y{n.,g)] _ ^^^^ 



+ ic.v^^^iF,(a + -;-;--^)j,(13) 

where is the confluent (Kunmier) hypergeometric 

function. 

The PDFs of Xo,fe, Y\ and Y2 are given by (12), when 
setting = Jls, a = 1 for the X^^,k\ = fii, a = MN for 
Yi; and O = Afili, a = N for the Yi, respectively. Thus, 
from (13), CFs of X, Y\ and Y^ are given by <l>x(w) = 
[$(cj,Og,l)]^^ ^yA^) = ^uj,ni,MN) and ^^y^M = 
^{LO,Mni,N), respectively. 

C. Outage Probability 

The CDF of the envelopes ratio G is determined as [9, Eq. 

(2)], 

FG{g)= / Fx{gy)fY{y)dy, (14) 
Jo 

where Fx{-) is the CDF of the desired signal envelope X. 
This CDF is expressible in terms of its CF by applying the 
Gil-Palaez theorem [10], 
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After introducing (15) into (14) and changing the orders of 
integration, we have 
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du> 



" ^x{oj)'i>Yi9'^) 
Im{$x(w)^r(fl'w)} 



doj , (16) 



where * denotes conjugate and Im{ } denotes the imaginary 
part of the argument. 

The straightforward approach is to estimate (16) by nu- 
merical integration. However, it is also possible to utilize 
an alternate approach, which will yield to an infinite series 
solution of (16). In [11] Beaulieu derived an infinite series for 
the PDF and the CDF of a sum of independent RVs, while 
[12] gives an alternative derivation that provided insights into 
the uses and limitations of the Beaulieu series. We use this 
alternative form of the BeauUeu series [12, Eq. (4b)], and 
express the CDF of X as 



Fxigx) = 1- 



2 lm{^xinuJo) exp(-jnwo5a;)} 
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The three alternative solutions for the system's outage proba- 
bility (i.e. the probability of SIR to fall below a given threshold 
z) are obtained by setting g = ^/z into (14), (16) and (18), 
which gives 
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Fx{y\/z)fY{y)dy , 
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d(j , (19b) 
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X $y(-(2n-l)a;oV^)} , (19c) 
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respectively. 

The exact OP can be derived for M = 2. The CDF of a 

sum of the envelopes of two Rayleigh-faded desired branch 
signals, Xo,i and Xo^2, is known [13], 
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where erf (•) is the error function. The derivation of the closed- 
form solution of (14) for the dual diversity case is provided in 
Appendix A, from which the outage probability at threshold 
z is determined to be 



Fz{z) = 1 



1 



1 + z/p 
X B 



[l + z/(2/3)]"+i/2 
1/2 1 1 
l + /3/^'2'"+ 2 



(21) 



where B( ;-, ) is the incomplete Beta function, defined by 

B^(a,6) = B{z;a,b) = f^'^l ~ t)'>-'^dt [14]. In (21), 
(a, P) = (2N, 7) for for incoherent interference combining, 
while {a, (3) = (N,'y/2) coherent interference combining, 
where 7 = fis /i^i represents the ratio of the average powers 
of the desired signal and a single interference signal in each 
diversity branch (also denoted as the average SIR per interferer 
per branch). 

In absence of diversity, M = 1, both interference combining 
scenarios converge and it is possible to directly solve (19a), 
which yields to the classic result for the OP, 



Fziz) = 1 
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(1 + ^/7)^ 



(22) 
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D. Average LCR 

In order to determine the average LCR of the random 
process G{t) by using (8)-(9), one typically needs to establish 
the joint PDF of the random processes G{t) and G{t) at any 
given moment t, ffjQ{g,g), as according to (10). However, 
we utilize an alternative approach, which circumvents explicit 
determination of fQQ{g,g). From (8)-(9), the time derivative 
of the envelopes ratio G is written as 



6 = lx-^y = ix-gi.. 



(23) 



Conditioned onY = y, the joint PDF fQQ{g,g) is calcu- 
lated as 

POO 

fGGi9,9)= fGG\Yi9,9\y)fYiy)dy, (24) 

where fy (y) is the PDF of the equivalent interference signal 
envelope Y. In (24), /(jgriyCff, ffly) is the conditional joint 
PDF of G and G given some specified value of the interference 
signal envelope Y ~ y, which is expressed as 

foGiYia^aly) = fciGyialg, y) ■ fcivigly) , (25) 

where fG\Y{9\y) is the conditional PDF of G given Y = y. 
Because of (8)-(9), it follows fG\Y{fj\y) = y ■ fx{gy), where 
fx{') is the PDF of the desired output signal envelope X. 

In (25), y) is the conditional PDF of G given 

some specified values of the envelopes ratio G = g and 
the interference signal envelope Y = y. Considering (23), 
this conditional PDF is determined as follows: Conditioned 
on G = g and Y = y, G is a linear combination of two 
independent RVs - the RV representing the time derivative of 
the desired signal envelope X{t) and the RV representing the 
time derivative of the equivalent interference signal envelope 
Yit). 

Under certain mathematical conditions, the envelope of 
the desired signal Xo,fc and its respective time derivative 
Xo.fe are independent RVs, and, at any given moment t, 
are characterized by the Rayleigh PDF and the zero-mean 
Gaussian PDF, respectively [1]. We conclude that the envelope 
of the desired signal at the EGC output X{t) and its time 
derivative X{t) are independent, since deriving (7) we get 



M 



(26) 
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Hence, X{t) is a zero-mean Gaussian RV with variance 
equal to the sum of the variances of the IID Gaussian 
RVs Xo^k{t) presumed to have equal powers, thus cr^ = 
M(T^ = (tt f „iqV Mils- This variance is valid for contin- 
uous wave (CW) transmission and two-dimensional isotropic 
scattering as according to the Clarke's model [1]. 

From (2) and (5), it is obvious that the instantaneous 
interference powers Y^ and Y2 can equivalently be repre- 
sented as the sum of MN and N IID squared Rayleigh 
RVs Rj with average powers f2i and Mill, respectively, as 
= ^^^^^^^ Rj. Finding the time derivative of both 
sides of the latter expression and specifying the values of 
the constituent Rayleigh RVs (thus fixing the values of Yi 



and Y2), one can easily conclude that both Yi{t) and 1^2 (i) 
are zero-mean Gaussian RVs with variances equal to cr^ and 
, respectively, independent of (t) and Y2{t) [15, Section 
3.i.l]. The latter conclusion is valid only if the variances of 
the time derivative of all constituent Rayleigh RVs are equal. 
In this case, assuming two-dimensional isotropic scattering, 
= {-KfrnifQ-i and <j\^ = {TrfmifMQi. 

Consequently, G is a zero-mean Gaussian RV with variance 
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^G\GY y2 '^X ^ yl'^Y ■ 



(27) 



Introducing (24) and (25) into (10), and changing the orders 
of integration, we obtain 

/•OO /"OO 

^G(5)=y^ 9'^9 fG\GYi9\9,y) fG\Yi9\y) fYiy)dy 

= fG\Y{g\y) fY{y)dy g fQ^Qy{g\g,y)dg . {22,) 
The inner integral in (28) is calculated by using (27), i.e., 



9fG\GYi9\9,y)dg = 



a 



G\GY 




Substituting (29) into (28) and considering /g|v'(!?|j/) = 
y ■ fxigy), we arrive at the important result for the average 
LCR of the envelopes ratio G at threshold g, 



Ng{9) = 




{gy)fY{y)dy. (30) 



The average LCR of G can also be evaluated in terms of 
the CFs of X and Y . Namely, after applying the Parseval's 
theorem over (30), we directly obtain 
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^x\ — \ ^Y {^) dio 



Re{$xH ^*Y{goj)}div , m 



where Rc{ } denotes the real part of the argument. 

The straightforward approach to obtain the average LCR of 
G is the numerical integration of (31). Alternatively, it is also 
possible to calculate the average LCR by using the infinite 
series solution after applying the Beaulieu series, similarly to 
the derivation of the OP. Namely, the PDF of X is expressed 
as [12, Eq. (4a)], 



fx{gy) = j; ^ Re{^x{nijJo)e^p{-jnijJogy)} + A2 , 

n=l,n odd 

(32) 

where A2 is an error term that tends to zero for large T, as 
assumed. Introducing (32) over (30) and changing the orders 
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of summation and integration, we obtain 

Ng{9) 



X 



27r 
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X ^^{nuog)} . (33) 



Thus, the three alternative solutions for the average LCR of 
the SIR Z at threshold z are obtained by setting g = into 
(30), (31) and (33), yielding 



Nz{z) 

Nz{z) -- 
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(34a) 
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Re{4>jc(w)$r(-WA/z)}(ia; , 
(34b) 
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X $r(-(2n-l)a;o\/i)}, (34c) 



respectively. 

An exact result can be obtained for the average LCR when 
M = 2. The PDF of the sum of the envelopes of two Rayleigh- 
faded desired branch signals, Xq^i and Xo,2, is known [13], 



X ( X \ I TT X ( X \ 

MX) = -exp ( -- ] + J— -exp ( -— j 



X erf 
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(35) 



The derivation of the closed-form solution of (30) for the 
dual diversity case is provided in Appendix A, from which 
the average LCR at threshold z for two-dimensional isotropic 
scattering is written as 
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1 4 -t- zal p~r(a + l/2) 
2^r 



1 



V 



Qs Via) 1 + 2/(2/3) 



1 («-l/2)-/J-l 



(1 -h 2//3)«-V2 



X B 



2 (l + z/(2/3))" 
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■ (36) 



l + l3/z' 2' 

In (36), (a,/3) = {2N,j) for the incoherent interference 
combining, and {a, (3) = {N,'y/2) for the coherent interfer- 
ence combining. 

If the desired and all interference signals are assumed to 
have same maximal Doppler frequency shifts, fmo = /mi. 
(36) is simplified into 

-r(a+l/2) 1 



Nz{z) = /mO\/7r- 



'z7/3 , [T (a - l/2)z//3- 1 




• (37) 



If the diversity is not employed at the receiver, M = 1, both 
interference combining scenarios converge, and it is possible 
to directly solve (34a) yielding to a well-known result for the 
average LCR at threshold z in Rayleigh fading when f^o = 
fmi [16, Eq. (17)], 



Nz{z) = /™oV27r 



T{N + 1/2) 



(38) 



T{N) (l + z/7)^ 
The AFD is calculated as Tziz) = Fz{z)/Nz{z). 

IV. Numeric Examples 

The numeric examples can be calculated either by numerical 
integration of (19b) and (34b) or by evaluation of the series 
(19c) and (34c). In this Section, we compare these two 
approaches/solutions in estimating the OP, the average LCR 
and the AFD, and then present some illustrative graphs. It is 
assumed that the desired and all interference signal have same 
maximal Doppler frequency shifts, fmO = fmi- 

The first approach requires utilization of a suitable numeri- 
cal integration technique embedded in the available computing 
software packages, such as MATHEMATICA. The adaptive 
Gauss-Kronrod quadrature (GKQ) method is particularly ef- 
ficient for oscillating integrands, such as those appearing in 
(19b) and (34b), which evaluates the integral at non-equally 
spaced points (abscissas) over the integration interval [17]- 
[18]. The GKQ method recursively subdivides the integration 
interval, reusing the abscissas from the previous iteration as 
part of the new set of optimal points, until the result converges 
to the prescribed accuracy. The number of abscissas (or, 
equivalently, the total number of integrand evaluations) and the 
number of recursive subdivisions needed to achieve the desired 
accuracy are not known ahead of computation, although the 
calculations of the new abscissas in each iteration based on 
[18] introduce a very low computational load. 

The second approach calculates the numeric examples by 
truncating the two infinite series solutions (19c) and (34c) to 
L non-zero terms. From the specialized solutions for the OP 
and the average LCR for M = 1 and M = 2, (21), (22), (37) 
and (38), one can conclude by induction that z and /3 appear 
together in the general expression for an arbitrary M, fornung 
the ratio z/(3. Consequently, the OP is calculated using (19c) 
and (13) as 

= ^ - ^ E imiimn - ivo, 1, i)r 

X $(-(2n - l)^'o 1, ")} . (39) 



n=l 



the normalized average LCR is calculated using (34c) and (13) 
as 



^ = ^ ^M + z/P Re{[$((2n - l)a;o, 1, 1)] 



M 



X ^-{2n-^)tJo^/z/0,^,a)}, (40) 

and the normalized AFD, ,fmX)Fz/Nz is calculated from the 
ratio of (39) and (40). Depending on the presumed scenario, 
{a, (3) = {MN, 7) - for the incoherent interference combining, 
and {a, (3) = {N,'y/M) - for the coherent interference 
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Fig. 1. First-order and second-order EGC output signal statistics vs. SIR Fig- 2. First-order and second-order EGC output signal statistics vs. SIR 
thresholds for various numbers of interferers when M = 3 thresholds for various diversity orders when N = 5 
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(b) Behavior of AFD 



Fig. 3. Second-order EGC output signal statistics vs. diversity order for 
various normalized SIR thresholds when N = 5 



Fig. 4. Second-order EGC output signal statistics vs. niunber of interferers 
for various normalized SIR thresholds when M = 2 



combining, while the normalized SIR threshold (NSIRth) is 
determined as 7/2 = Os/(^i2^)- Note that (39) and (40) actu- 
ally estimate (19b) and (34b), respectively, by sampling their 
integrands at equally spaced abscissas, whereas their number 
L and locations (odd multiples of loq) are given ahead of 
computation. There is a tradeoff between the absolute accuracy 
and the selection of T and L. A larger value of T results in 
greater accuracy, but more nonzero terms L must be used [11]. 

Using MATHEMATICA, we compared the computational 
burden between the two approaches/solutions by calculating 
the same numerical examples with same prescribed absolute 
accuracy of ±10~*. In utilizing the first approach, we set 
the target accuracy into the computing software, and, for a 
given set of input parameters (NSIRth, M and N), obtain 
the integration result and the respective value of the built-in 
variable that counts the number of integrand evaluations. In 
utilizing the second approach, the values of T and L needed to 
achieve the target accuracy are obtained empirically by trying 



multiple combinations of T and L, thus yielding to typical 
values of T between 40 and 100, and L - between 100 and 
200. 

For a given selection of T and L, (39) has a significantly 
better rate of convergence then (40), so the accuracy of (40) 
determines the accuracy of both the normalized AFD and 
normalized average LCR. Similarly, the GKQ of the OP (19b) 
produces a numerical result with better convergence rate and 
less integrand evaluations compared to the GKQ of the average 
LCR (34b). 

We estabUshed that the first approach introduces higher 
computational load, thus requiring longer computation times. 
Namely, for the range of the input parameters shown on 
Figs. 1, 2 and 3, a single numeric integration using the GKQ 
method requires between 100 and 500 integrand evaluations 
to estimate a single value of the OP or the average LCR. In 
the same range of the input parameters, the truncated BeauUeu 
series typically require fewer number of integrand evaluations 
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(L nonzero terms) to achieve the same accuracy, thus yielding 
to shorter computation times to obtain the respective results. 

It is also observed that the increase of a adds to the com- 
putational burden of the GKQ of both (19b) and (34b), since 
their integrands become more rapidly oscillatory (particularly 
emphasized for the average LCR calculations), thus requiring 
more integrand evaluations (increasing toward 500) to achieve 
the desired accuracy. The number of required nonzero terms 
i in the BeauUeu series increases (toward 200) for lower 
NSIRth (i.e. higher threshold z), but rapidly decreases with 
the increase of the diversity order M. Compared to the first 
approach, however, the computational loads of the Beaulieu 
series in estimating OP and average LCR are less dependent 
from their input parameters. 

Figs. 1 and 2 depict the OP, the normalized average LCR 
and the normalized AFD versus the NSIRth, with N and M 
appearing as curve parameters (M = 3, = 1,5 and 10 in 
Fig. 1, and M = 1, 2 and 5, A'' = 5 in Fig. 2). Note that if 
NSIRth < dB then z > Os/^^i, while if NSIRth > dB 
then z < Jls/f^i- 

As expected. Figs, la and 2a show that the OP is a 
monotonically decreasing function from the NSIRth, and that 
its values match [2, Fig. 1] and [3, Figs. 2 and 3] in the 
respective NSIRth ranges for given N and M. Figs, lb and 
2b show that the average LCR reaches its maximum for some 
specific NSIRth, whose value depends on the selection of M 
and N. Figs. Ic and 2c show that the AFD decreases by the 
increase of the NSIRth (i.e. by the decrease of the threshold 
z). 

It is also obvious that, for a given values of NSIRth, M 
and A'^, the average LCR and the AFD curves for incoherent 
and coherent interference combining almost coincide when the 
NSIRth is above the value that maxinuzes the average LCR, 
while they differ below this value. 

Figs. 3 and 4 depict the influence of the diversity order M 
and the number of interferers N over the average LCR and the 
AFD for different values of NSIRth. Depending on whether 
NSIRth is set below or above the value that maximizes the 
average LCR, the average LCR may increase and/or decrease 
by increasing M (Fig. 3a), while the AFD monotonically 
decreases (Fig. 3b). The average LCR may also increase 
and/or decrease by increasing N (Fig. 4a), while the AFD 
monotonically increases (Fig. 4b). 

Extensive Monte Carlo simulations conducted in MATLAB 
have validated all numeric examples presented in this Section. 



diversity branch. The infinite series solutions were determined 
for an arbitrary diversity order after successively applying the 
CF method, the Parseval's theorem and the Beaulieu series 
over the integral expressions for the OP and average LCR. 
Compared to the numerical integration method commonly im- 
plemented in the computing software packages, we concluded 
that the Beaulieu series solutions introduce less computational 
burden and minor sensitivity to the input parameters. For 
the dual diversity case, the average LCR and the AFD were 
determined as exact closed-form solutions in terms of the 
gamma and the beta functions. 

For the interference-Umited EGC systems, the desired 
branch signals coherently combine, while the branch signals 
from each interferer can combine either coherently or inco- 
herently. Our analytical solutions incorporate both combining 
scenarios, yielding to somewhat different numeric values for 
the average LCR and the AFD. The differences are more 
evident when the SIR threshold is set above the average SIR 
per interferer per branch. 

One can further alleviate the assumption for the equal 
branch powers of the desired signal replicas. It is straightfor- 
wardly obvious from (26) that same analytical approach is also 
applicable for determination of the average LCR in the case 
of unequal branch powers of the desired signal. The derivation 
of the respective solutions is trivial and omitted in this work. 



Appendix A 

Introducing the expressions for the CDF of X (20) and the 
PDF of Y (12) into (14), and also using [14, Eq. 3.478(1), 
6.286(1)], one obtains 



Fcig) = 1 



1 



where 2-Fi(a; b; c; z) is the Gaussian hypergeometric function 
[14]. We then successively apply transformations [14] 



2Fi{a;b;c;z)= { -;^] 2F1 { c - a;b;c; ] (A.2) 



I- z 



z-1 



and 



V. Conclusion 

This paper derived the analytical expressions for the average 
LCR and the AFD of coherent EGC wireless communication 
systems subject to CCI and Rayleigh fading. 

The solutions for the average LCR were derived by a novel 
analytical approach that circumvents the necessity of finding 
the explicit expression for the joint PDF of the instantaneous 
SIR and its time derivative. They have been expressed in 
forms of an infinite integral solution and an infinite series 
solution, assuming IID equal-powered interference signals' 
replicas and IID equal-powered desired signal replicas in each 



2Fi{l;b;c;z) = ^{l-zy-^-' B{z;c-l,b-c+l) (A.3) 



over (A.l) and obtain 



Faig) = 1 



1 + 52//?; [l + .g2/(2/3)]"+i/2 



The result (21) is obtained directly from (A.4) by setting g = 



9 



Introducing the expressions for the PDF of X (35) and the 
PDF of Y (12) into (30), and also using [14, Eq. 3.478(1), 
6.286(1)], one obtains 



Ng{9) = 



1 r(a + l/2) 



2ti y\l T{a) 
g^/^ , r^r(a + 3/2) 



(l + <7V/3)"+^/^ V r(a) 

^ (g/V?)^ p /I 3 3 1 
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(1 +.g2/(2/3))"+3/2 ' 'V2'2 '2' 1 + 2^/g2 

rr(« + i/2) 



r!s r(a) (l + g2/(2/3))"+i/2 

/^l 1 3 1 
X 2i^i 77; 77 + a; -; 



2' 2 '2' 1 + 2/3/52 



(A.5) 



After successively applying both (A.2) and (A. 3) over (A.5), 
one obtains 



r r^r(a + i/2) 
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1 r(a + 3/2) 
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,l + /3/g2'2 
We further apply the following identity [14] 



z°'(l~z)'' b 
Biz;a,b+1) = ^ + -B{z;a,b) (A.7) 



a + b a + b 



over (A.6), and obtain 
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The result (36) is obtained directly from (A.8) after setting 
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